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Echoes of a decaying planetary system: the gaseous and dusty 
disks surrounding three white dwarfs 

C. Melis^'2, M. Jura^ L. Albert^ B. Klein^ B. Zuckerman^ 
cmelisOucsd . edu 

ABSTRACT 



We have performed a comprehensive ground-based observational program 
aimed at characterizing the circumstellar material orbiting three single white 
Q> ■ dwarf stars previously known to possess gaseous disks. Near-infrared imaging 

C/^ ■ unambiguously detects excess infrared emission towards Ton 345 and allows us 

1^ ■ to refine models for the circumstellar dust around all three white dwarf stars. 

We find that each white dwarf hosts gaseous and dusty disks that are roughly 



jh ■ spatially coincident, a result that is consistent with a scenario in which dusty 

^ . and gaseous material has its origin in remnant parent bodies of the white dwarfs' 

planetary systems. We briefly describe a new model for the gas disk heating 

mechanism in which the gaseous material behaves like a "Z 11" region. In this 

Z II region, gas primarily composed of metals is photoionized by ultraviolet light 

CN . and cools through optically thick allowed Ca Il-line emission. 

O 

(N 

[■^ ■ Subject headings: circumstellar matter — planet-star interactions — stars: in- 

dividual (Ton 345, SDSS J122859.93+104032.9, SDSSJ104341.53-f085558.2) — 
white dwarfs 



C^ ■ 1. Introduction 

White dwarfs represent the endstate of stellar evolution for the majority of stars within 
our galaxy. Most stars with mass less than ~8 Mq will end up as these exposed cores of 
nuclear ash that slowly cool over time as they radiate away their energy. The planetary 
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systems that may have orbited these white dwarfs while they were on the main sequence 
have been shown to be able to suryive through the asymptotic giant branch (AGB) phase of 



stellar evolution (jSackmann et a/.lll993l : iDuncan fc Lissauerlll998l : iDebes &: Sigurdssonll2002 



Burleigh et a/.ll2002l ). What is the ultimate fate of these remnant planetary systems? 



The work of iDuncan fc Lissauen ( 1l998l ) showed numerically that giant planets orbiting 
white dwarfs descendent from massive progenitor stars are likely to become dynamically 
unstable on timescales less than ~10 Gyr. Less massive progenitor stars are more likely to 
have giant planets with stable orbits at the end of their stellar evolution. However, these less 
massive progenitor stars are also expected to obliterate any terrestrial-like planets orbiting 
within ~1 AU during their AGB evolutionary phase, wherea s all terrestrial planet s around 



more massive stars will likely survive through the AGB (e.g., iSackmann et a/.lll993l ). Rocky 
objects; hke comets, asteroids, or terrestrial planets; are expected to survive to the white 
dwar f evolutionary phase for a variety of initial sizes and orbital semi-major axes (jjura 
20081 ). but may suffer orbital de-stabili zation through the perturb ations of giant planets and 
be flung towards the white dwarf star ( JDebes fc SigurdssonI 120021 ). 



Once within the white dwarf's Roche radius, rocky bodies will be tidally shredded 
and eventually accrete onto the white dwarf's photosphere ([Debes fc SigurdssonI l2002l : IJura 
2003al : IZuckerman et a/.l 120071 : iKlein et a/.ll2010l : iDufour et al\\201(j) . Sub-micron size grains 
orbiting within the Roche radius of their host white dwarf star can produce a noticeable 
infrared excess above what would be expected from the white dwarf photospheric emission 
alone. Directly detecting rocky objects in orbit around white dwarfs will be exceedingly 
difficult; in contrast, the dusty finale of rocky objects that fall into a white dwarf's Roche 
radius and the eventual atmospheric pollution from the residual rocky object material leave 
a tell-tale sign that there is a dissolving remnant planetary system. Thus, much like current 
searc hes around main sequence stars utilize debris disks as signposts for planetary systems 



(e.g., IZuckerman fc Songil2004J : iKrivovl 120101 ). remnant rocky planets around white dwarfs 
are currently best investigated indirectly through searches for white dwarf debris disks and 
atmospheric pollution. 

Past searches for remnant planetary systems around white dwarfs have focused on direct 
and indirect detection of giant planet companions (e.g., iHogan et a/.ll2009t iMuUally et al. 



20081 . and references therein) and the detection of dust contained within a white dwarf's 



Roche radius as discussed above (e.g., iFarihi et a/.ll2009l and references therein). To date 



searches for giant planets have had little success (but see iMuUally et al\ l2008l . 120091 ) . The 
first evidence for a remnant planetary system around a white dwarf came from the work of 
Zuckerman &: BecklinI (119871 ) in which it was discovered that the ne arby white dwarf star 
G29— 38 hosted infrared excess. This infrared excess was interpreted by lZuckerman &: Becklin 
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( 119871 ) to probably be the signature of a brown dwarf companion. However, subsequent ob- 
servations of G29— 38 found evidence for a silicate emission feature at ~10 /xm (JGraham et al. 



1990l : iTokunaga et a/.lll990l : iReach et a/.ll2005l ) indicating that the excess originated from a 
dusty debris disk. Furth er ground-based investigation of whit e dwarfs decades later reveale d 
three new debris disks JSecklin et a/.lboosi : Ixihc et all^OO^ . bood : iKihc fc Redfieldl boOTh 



It to ok the launching of the mid-infrared optimized Spitzer Space Telescope (IWerner et al. 



2004 ) to obtain the necessary sensitivities to probe for a statistical sample of debris disk- 
hosting white dwarfs. Spitzer proved excellent at discoverir ig debris disks arou nd white dwarf 



stars, increasing the sample size by nine (see Table 1 in iFarihi et all |2009| . and references 
therein). Spitzer observations that reached white dwarf phot ospheric flux l evels enabled a 
statistical measure of the incidence of white dwarf debris disks. iFarihi et al\ ( 120091 ) compiled 
all previous Spitzer observations and found that ~l-3% of all single white dwarfs with cool- 
ing ages <0.5 Gyr (Teg^lO, 000-20, 000 K) have debris disks and that ~50% of white dwarfs 
with implied metal accretion rates of M >3xl0^ g s~^ have dusty material orbiting them. 
All white dwarfs discussed above have debris disks contained almost entirely within their 
respective white dwarf's Roche radius indicating disrupted rocky objects. 



Contemporaneous with the Spitzer white dwarf disk searches, iGansicke et al\ (120061 ) 
noted pecuhar double-peaked gas emission lines from a ho t (Tpff=22,000 K ) white dwarf 
observed as part of the Sloan Digital Sky Survey (SDSS; lYork et al\ |2000[ ). This single 
white dwarf, dubbed SDSS 1228-1-1040 (spectral type DAZ), was interpreted to have an 
orbiting gaseous metal disk that had its origin in disrupted asteroidal material. Searches 
through the rest of the SDSS white dwarf spectroscopic database unc overed two more gas- 
disk w hite dwarfs (SDSS1043+0855, DAZ, and SDSS0843+2257, DBZ: lGansicke et alhm 
l2008al . hereafter SDSS1043 and Ton 345). Spitzer and ground-based observations of SDSS 
1228-|-104(y (hereafter SDSS1228) show that this gas disk-hosting white d warf is also orbited 



by du sty material that coincides positionally with the gaseous material (JBrinkworth et al. 



20091 ) . Such a discovery heralds the use of such gas and dust disk hosting white dwarfs as 
benchmark objects for studying the decay of planetary systems as the two disk components 
yield complementary insights; for example, gas disk measurements can directly inform one 
about the dynamics of the disk while dust disk modeling indirectly provides information 
about disk structure. 

Through judicious observing with both ground- and space-based facilities there is now 
convincing evidence that white dwarfs show signs of dissolving planetary systems. However, 



^ Farihi et ali (|2010f ) performed AKARI observations of Ton 345 that indicate there is excess infrared 
emission at 2.3, 3.2, and 4.3 /ini; these observations are contaminated by a nearby background object 
discussed in Section [2l 
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the sample size of debris disk white dwarfs is still small and would benefit from additional 
objects. Furthermore, it is not yet clear how the gas disk white dwarfs fit into the dissolving 
planetary system framework that explains the white dwarf debris disk observations. Thus, in 
an effort to better understand the final throes of planetary system evolution, ground-based 
ne ar-infrared and optical spectroscopi c observations of the three gas-disk systems discovered 
by iGansicke et al\ (120061 . 120071 . l2008al ) were performed and are reported herein. 



2. Observations 

2.1. WIRCam Imaging at the CFHT 

Initial near- infrared (NIR) images of Ton 345 and SDSS 1043 were obtained with WIRCam 
mounted on the 3.6-m CFHT telescope (jPuget et a/.ll2004f ) on the nights of UT 24 February 
2008 and UT 19 June 2008, respectively. We used a 6- to 9-point dither pattern with 60" 
offsets between dithers. For Ton 345 we took exposures in the J, H, and Kg bands with 45, 
15, and 25 seconds per dither (yielding total integration times of 270, 630, and 450 seconds, 
respectively). The H-band observations for Ton 345 were repeated on UT 25 February 2008 
in better sky brightness conditions. Both H-band data sets yielded the same result. For 
SDSS1043 we took exposures only in the Kg-band; these observations were performed in 
a manner similar to that described for Ton 345. Images were preprocessed, non-linearity 
corrected, and sky subtracted at CFHT with the T'iwi pipe li nci a nd media n-stacked usin g 
the suite of Terapix software (sextractor - iBertin k. Arnoutd (119961 ). scamp - iBertinI ( 120061 ). 
and swarplj) with the 2MASS point-source catalogue as the reference for external astrome- 
try. Photometry of the two white dwarfs and all 2MASS sources was performed using the 
FLUX_AUTO measurement in SExtractor, version 2.4.4. Absolute photometry was anchored 
on the 2MASS system using a weighted average of all (>200) 2MASS point sources in the 
full field of view and 5 iterations of 3cr outlier rejection. The uncertainty on this calibration 
is ~0.03 magnitudes in all three bands for objects with fluxes similar to Ton 345 and ~0.05 
magnitudes in Kg for objects similar in flux to SDSS1043. 

WIRCam data for Ton 345 are summarized in Table [T] which also includes SDSS and 
GALEX fluxes. Images of Ton 345 and a ~1.8" separation background source, presumably a 
galaxy, are displayed in Figure [H Large beam instruments (like Spitzer or AKARI) would be 
contaminated by this red background object. The WIRCam result for SDSS1043 is presented 



^http: //www. cfht .hawaii . edu/Instruments/Imaging/WIRCam/IiwiVersionlDoc .html 
''http: //www. astromatic .net/ 
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in Table H 



2.2. NIRI Imaging at Gemini-North 



Follow-up observations of Ton 345 were obtained with NIRI ( iHodapp et a/.ll2003l ) mounted 



on the Gemini North Telescope. Observations were performed on UT 24 May 2008 during 
Director's Discretionary Time (program GN-2008A-DD-5). 5-point dither pattern image sets 
were repeated for J, H, and Kg with total integration times of 105, 120, and 150 seconds, 
respectively. A 4-point dither pattern was repeated for L' until we accrued a total integration 
time of 2520 seconds. 

After pre-processing to correct for the known NIRI non-linearitjQ, data were reduced 
using in-house IDL software routines. For each filter, science frames were median combined 
to generate a "sky-median" frame which was then subtracted from each science frame. Sky- 
subtracted frames were then fiat-fielded using exposures of the illuminated telescope dome 
for J, H, and Kg. L' frames were flat-fielded using the median stack of all science frames 
normalized to unity. Reduced science frames were registered by reversing the header-recorded 
values for the telescope offsets between dithers. 

Photometric measurements for Ton 345 were calibrated by comparing extracted fluxes 
for each filter to standard stars fluxes (Leggett et al 2003, 2006) observed and reduced in a 
similar fashion. For JHKg the standard star was FS 15, for L' it was HD 84800. We note 
that HD 84800 appears to have a ~0.4" extent in our L' images; no telescope problems were 
noted at the time of observation. JHKs fluxes were extracted with an aperture that yielded 
~85% encircled energy (with a negligible difference between Ton 345 and FS 15). For L' flux 
calibration we used a 2.0" diameter aperture to extract the flux of HD 84800 and Ton 345. 
Uncertainties for JHKg fluxes were determined by examining the dispersion of the extracted 
fluxes for each of the individual reduced science frames. Since this was not possible in our L' 
imaging sets, we instead made ten randomly placed extractions on the final, median-stacked 
L' image using the same aperture that was used to extract the flux for Ton 345. The standard 
deviation of this set of values was adopted as the flux measurement uncertainty. The final 
signal-to-noise ratio (S/N) is ~50 for JHKs while for L' it is ~3-5. 

NIRI photometric results are listed in Tabled! Figures [T] and H] show the NIRI images of 
Ton 345 and the ~1.8" separated galaxy. Table [3] lists fluxes for the galaxy extracted using 
similar apertures and procedures as for Ton 345. 



"'http : //staff . gemini . edu/~astephens/niri/nirlin/ 
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2.3. Gemini Imaging at the Shane 3-m 



Observations of SDSS1043 in the J-band were perform ed UT 10 October 2008 with tlie 



Gemini Twin- Arrays Infrared Camera (JMcLean et a/.lll993l ) mounted on the 3-m Shane tele- 



scope at Lick Observatory. We used a 5 position dither pattern with 20 second integrations 
of 5 coadds per position resulting in 500 seconds total on source integration time. The ~3' 
field-of-view of the Gemini instrument enabled simultaneous observations of two 2MASS 
stars for use in flux calibration. Although Gemini hosts two infrared cameras, data for the 
longer wavelength chip were unusable due to instrumental difficulties. 

Data were reduced and fluxes extracted much like that described above for the Ton 345 
NIRI data. Exceptions include no non-linearity correction, the use of twilight fiats instead 
of dome fiats, and image registration using bright point sources within the field. The Gemini 
J-band fiux for SDSS1043 is reported in Table [2l 



2.4. HIRES Spectroscopy 



Keck HIRES (jVogt et a/.lll994l ) optical echelle spectra were obtained for each of the gas- 
disk white dwarfs. Table H] lists observation dates and HIRES instrumental setups. Data were 
reduced via two software packages: IRAF and MAKEE. Standard IRAF tasks were used in 
the reduction and extraction of the HIRES data. After reduction and extraction, low order 
polynomials were fit to each order to bring overlapping order segments into agreement. 

The output of each software package was compared to ensure no reduction artifacts 
were contaminating gaseous emission fines. We found that additional processing steps were 
necessary for data taken in the red setup. In some observing setups blue light in second order 
leaked into the ~8200 to 9000 A range. This could potentially "veil" absorption and emission 
lines in this region, especially the Ca II A8662 transition. To estimate the contaminating 
flux present we extrapolated the flux from blue contaminant orders that fell in between red 
orders. We believe that this excess fiux was removed for most observations, only the 2008 
November observations of SDSS1228 and SDSS1043 might have significant residuals due to 
strong Balmer absorption lines that were difficult to model. 
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3. Results 

3.1. Imaging 

Table [T] lists measured near-infrared fluxes for Ton 345 from the CFHT and Gemini 
Nort h data sets. To place these near-infrared measurements in context we queried SDSS 



DR7 ( lAbazajian et a/.ll2009l ) and GALEX for ugriz and ultraviolet fluxes, respectively. From 
the measured magnitudes alone it is apparent that Ton 345 hosts signiflcant near-infrared 
excess (see Section H]). 

Table [2] lists GALEX ultraviolet fluxes, SDSS ugriz, and the ground-based near-infrared 
photometry for SDSS1043. Our data alone do not conclusively show that there is infrared 
excess emission towards SDSS1043. But, when these ground-based near-infrared measure- 
ments are combined with Spitzer photometry at longer wavelengths (C. Brinkworth et al. 
2010, in preparation), signiflcant excess infrared emission is apparent (see Section |4]). 



3.2. Optical Spectroscopy 

Double-peaked metallic gas emission lines are detected towards all three white dwarfs. 
Emission lines from Ca 11 and Fe 11 are shown in Figures [3] and HI Ca 11 H & K emission 
lines are detected only for SDSS1228 (Figure |3] and Table [5]). Fe 11 A5169 emission lines are 



detected for SDSS1228 and Ton 345 (F igure H Figure 



lines for SDSS1228 that were noted by lGansicke et al\ (|2006[ ) are not detected. Measurable 



and Table [5]) . Fe 11 A5018 emission 



HI and 
1 12006 



quantities from emission features in the spectra, and associated derived quantities, can be 
found in Tables E] and IB We defer quantitative analysis of absorption feature equivalent 
widths and elemental abundances to future publications. 

Table [5] contains information on the velocity separation between emission peaks of the 
same transition, the full velocity width at zero power of each double-peaked emission feature, 
the "peak midpoint" , and the maximum velocity gas detected in the blue and red wings of the 
double-peaked emission features (vmaxSim). We estimate the wavelength centroid position 
of emission peaks by fltting either Gaussian or Lorentzian proflles over a ~10-15 A range 
centered on the emission line peak. Three to four flts per peak are performed with the 
continuum anchor points for the fltting function at different locations within the noise of 
the spectrum. The peak centroid is taken to be the average of these multiple flts while an 
estimate of the uncertainty on the centroid comes from the standard deviation of these flts. 
The peak midpoint position is then determined by averaging the two emission peak centroids 
for each transition. The peak centroid and peak midpoint wavelength values are converted 



into velocities relative to the rest wavelength for the corresponding transition they are being 
measured from. 

Measured radial velocities for each star are reported in the second column of Table O 
We compute Galactic UVW space motions using the gravitational redshift-corrected radial 
ocities (where gravitational redsh ifts in Table |6] are computed using the stellar parameters 



ve 
of 



Gansicke et al 



20061 . 120071 . l2008al . which are reproduced in Table [6]), 2MA SS skv positions 



(ICutri et a/.ll2003l ). proper motion information from the NOMAD catalog (jZacharias et al. 



20041 ). and distances computed from the Rwp/ Dist parameter f rom our model fits (see Section 



Sland Table [7]) and the gravity values given in lGansicke et al\ ( 120061 . 120071 . l2008al ) . We adopt 
an uncertainty of 10% for the distances. 



3.2.1. Gas Disk Evolution 

With m ultiple epochs o f spec t ra, and by c omparison to the Ca II IRT emission lines 
presented in iGansicke et al\ (120061 120071 . |2008aJ) , we can probe variability in emission line 
strength and morphology. Table [5] contains emission line equivalent widths, velocity sep- 
aration between emission peaks of a given transition, and the full width at zero power of 
emission lines for each epoch a gas emission line was observed. Figure S] shows multiple 
epochs of Fe II A5169 emission lines for SDSS1228 and Ton 345 and multiple epochs of Ca II 
IRT emission in Ton 345. There are no significant variations detected in the Fe II emission 
lines (but see caption to Figure H]) . 

Ton 345's Ca II IRT emission line equivalent widths and morphologies appear to agree to 
within the errors ass ociated with each indiv idual HIRES epoch. This is in contrast to the two 
epochs presented in iGansicke et al\ ( l2008al ) in which the equivalent width of the sum of the 
three Ca II IRT lines decreases by a factor of ~1.5 and a significant morphological change in 
the emission profile is seen between 2004 Dec and 2008 Jan. The measureme nts for the sum 



of the Ca II IRT emission equivalent width reported herein agree with the IGansicke et al. 



( l2008a[ ) 2008 Jan value to within the la errors. Although speculative, this might suggest 
that an episodic disk feeding event occurred before the 2004 Dec epoch and that the 2008 
measurements are probing a "quiescent" phase. 



m 



Comparing the HIRES ob servations of SDSS1228's Ca II IRT complex to those presented 



Gansicke et al\ (120061 120071 ) suggests the stronger of the two emissi on peaks has switched 



from the red side of the double peaked emission complex (as seen in IGansicke et al\ 12006 
20071 ) to the blue side (HIRES data). This parity switch in the peak streng th is confirmed 



when we degrade the HIRES data to the resolution of the spectra presented in lGansicke et al. 
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( 120061 ) ■ Interpretation of this par ity switch is co n aplica t ed by the sparse time samphng of 
the emission hne complexes. The iGansicke et all (l2006l . 120071 ) spectra have epochs of 2003 
Mar (SDSS) and 2006 Jul (WHT). The morphological appearance of the Ca II IRT emission 
lines in the SDSS and WHT spectra suggests little change between these two epochs. It 
is not clear what would change the parity in the emission peak strength between the 2006 
WHT epoch and the 2008 HIRES epoch. One possibility is that the disk is clumpy and we 
have witnessed the clumps orbiting around SDSS1228. Additional spectroscopic monitoring 
of SDSS 1228 is necessary to confirm and characterize such orbiting clumps. 



4. Dust Disk Parameters 



Ton 345 's photometry suggests near-infrared emission in excess of what one would ex- 
pect from the photosphere of the white dwarf st ar alone. The comprehensive data set in 
Table d] is modeled as a flat, passive, opaque disk (|Jurall2003bl : I Jura et a/.ll2007al ) orbiting a 
~18,600 K effective temperature white dwarf. Ton 345's photospheric emission is assumed to 
be well represented by a blackbody curve; such a model is reasonable for Helium dominated 
atmosphere white dwarfs that do not exhibit strong ultraviolet flux suppression from Hy- 
drogen absorption. Figure [5] displays the white dwarf and disk model while Table [7] reports 
the model parameters. The photospheric and excess fluxes are reproduced succesfuUy by 
the white dwarf and disk model. However, because the measurements used for Ton 345 only 
extend out to wavelengths as long as L', the temperature at the outer edge of the dust disk 
is observationally unconstrained. In the model shown in Figure [5] and Table [7] we assume 
Touter of the disk is 1000 K. Touter could be, and likely is, smaller than this value (see be- 
low). Further observational constraint of the disk outer edge temperature will require longer 
wavelength data that is uncontaminated by the nearby background object. 

In the interest of comparing a homogeneous set of disk model parameters, we reflt 
the data for SDSS1228 presented in ISrinkworth et all J2009h and the data for SDSS1043 
presented in C Brinkworth et al. (2010, in preparation) in the same manner as Ton 345. 
These model flts and measured fluxes are shown in Figures [6] and [7] with model parameters 
reported in Table [71 It is noted that SDSS1043's outer disk temperature is poorly constrained 
even with the inclusion of IRAC photometry. In an effort to illustrate the range of viable 
model parameters that can flt the measured infrared excesses, we plot in Figure E] three 
models for SDSS1228 that have different dust disk parameters (Table [7]). The resultant 
parameter extremes reported in Table [7] are taken to be the ~5cr confldence limit on each 
model parameter. We then derive the la model parameter uncertainties by dividing the 
difference between parameter extremes by flve; these uncertainties are reported along with 
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SDSS1228's best fit model parameters in Table [71 



For the case of optically thick, fiat dust disks Eq. 1 from iJural (l2003bl ) can be used to 
determine at what radial distance from each white dwarf dust particles of certain temper- 
atures reside. As previously mentioned, for Ton 345 and SDSS1043 the temperature of the 
outer region of their orbiting dusty material is poorly constrained by observations. Rather 
than take the model outer disk temperatures for Ton 345 and SDSS1043 (which are better 
interpreted as upper limits), we instead estimate their outer disk temperatures under the 
assumption of an asteroidal debris model for the origin of the dusty dis k. Such a scenario re- 
quires a rocky obj ect to fall within the Roche radius of the white dwarf (IDebes fc Sigurdsson 
2OO2I : |jurall2003bf ). The following expression for the Roche radius is used ( DavidssonI Il999l ) : 



R. 



tide 



a 



tide 



pa 



1/3 



-R* 



where within Rude (the radial separation from a star of radius i?* and density p*) a rocky 
object of density pa would be disrupted. The factor Cude is a constant of order unity that 
is determined by the tensile s trength, shape, rotation r ate, and orbital parameters of the 
rocky object (lDavidssonlll999t iHolsapple fc Michelll2008l ). It is unlikely that all objects that 
were tidally shredded around SDSS1228, Ton 345, and SDSS1043 had the same physical 
parameters; however, in the absence of any such information it is assumed that this is the 
case. Such an assumption can potentially introduce a significant systematic offset in the 
derived outer dust disk radii. Using the well-modeld data for SDSS1228 (Table [7]) and 
assuming an asteroid density of ~3 g cm~^ a value of 1.45 for Cude is estima ted. Combining 



these values with th e mass and radius of the white dwarfs (see Table |6]) from lGansicke et al. 



( I2OO6I . I2OO7I l2008al ) we calculate the tidal radii for each white dwarf. These values for Ton 
345 and SDSS1043 are reported as ^outer,dust in Table El The value reported for SDSS1228 
in Table JHl is the disk outer radius as derived from fitting its Spitzer data. 



5. Gas Disk Parameters 



Using the HIRES optical spectra we can estimate the radial extent of the gaseous disks 
orbiting the three white dwarfs. This is done by measuring the highest velocity gas emission 
in the HIRES spectra (vmaxSini in Table JS]), emission that corre sponds to the innermo st orbit 
of the emitting gas-phase metals around the white dwarf (e.g.. lHorne fc Marshlll986[ ). How- 
ever, this value is degenerate with the inclination angle of the disk, z, where the degeneracy 
is fsini {i of 0° would correspond to a face-on disk). The dust disk model fits (Section H]) can 
help constrain the disk inclination angle assuming the dust and gas disks are coplanar. It is 
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noted that the inchnation angles as derived from the dust disk fitting can vary by ~20-30% 
from the best-fit value (5cr; see Tableland Section H]). Line-of-sight gas disk velocities for 
the three white dwarfs are corrected using the dust disk model best-fit inclination angles. 
Assuming Keplerian orbits of the gas and dust disks, we use the w hite dwarf masses and 
radii (see Table |6]) as modeled by iGansicke et all (120061 . 120071 . l2008al ) to derive the gas disk 
inner radius in units of white dwarf radii. We obtain the full diniensions of the g aseous disks 
by combining these measurements with the lGansicke et all (120061 . 120071 . l2008al Jbl) models (see 
Section E]) in which the gas disk outer radii have been estimated. These results are reported 
in Table |8] and illustrated in Figure [HI 

We note the discre pancy between ou r infer red gas disk inner radius for SDSS1228 and 
the same as quoted in iBrinkworth et all (120091 ) . The origin of this discrepancy has to do 
wi th two factors. The first factor is responsible for the smaller gas disk inner radius quoted 

^27 Rwn — and results from use of the maximum 



m 



Brinkworth et al\ (l2009f ) - R 



dinner, gas 



Ca II infrared triplet (IRT) gas emission velocity of ~1270 km s ^ as reported in Table 1 in 



Gansicke et al\ (120061 ). Re- measurement of the spectra presented in iGansicke et al\ (120061 ) 



provides a m aximum gas ernission velocity more in line with that measured herein (Table 
use of the IGansicke et al\ l2006l . erroneous VmaxSini will reduce the gas disk inner radius 
by a factor of ~1.8 relative t o our inferred gas dis k inner radius) and what is quoted in 
the supplementary material of IGansicke et al\ ( 20061 1^1. The second factor is a 20% increase 
in the inner disk radius relative to that presented herein (wher e the ratio of 1.8/1.2= 1.5 is 
the ratio between our inferred gas disk inner rad ius and that of IBrinkworth et a/.ll2009[) and 
comes from a matter of interpreting Figure 1 of iHorne &: MarshI (| 19861 ). The inner gas disk 
radius repor ted in Table [8]f o r SDS S1228 agrees with the value quoted in the supplementary 
material to iGansicke et al\ (2006[o. 



5.1. Stable Eccentric Disk at Ton 345 

Here a novel method of estimating the gravitational redshift for the gas disk-hosting 
white dwarf stars is attempted. This method employs the gas disk emission lines (from 
which we seek to measure the white dwarf systemic velocity) and photospheric absorption 
lines (which contain velocity components from the systemic velocity and the gravitational 
redshift from the compact white dwarf). To estimate the gravitational redshift we measured 
the photospheric absorption line radial velocity for each white dwarf (see RVobs in Table 



^The values reported in the main article text and tables of lGansicke et al\ ( 20061) are different from those 
reported in the supplementary online material. 
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E]) and the systemic radial velocity for each white dwarf-disk system from gaseous emission 
lines (see "Peak Midpoint Velocity" in Table |5] and discussion in Section [3^ . Our estimated 
gravitational redshift for SDSS1228 (RVobs — [Peak Midpoint Velocity]avg= 55±5 km s~^; 
based only on the Ca II IRT emission lines) reproduces well the expected gravitational 
redshift (44 km s~^, see Table [6]). Estimates for SDSS1043 are inconclusive due to the low 
S/N in the gas emission lines. However, our estimate for Ton 345 (RVobs — [Peak Midpoint 
Velocity]avg= 131±6 km s~^; based only on the Ca II IRT emission lines) is quite different 
from the expected gravitational redshift (44 km s~^). 

Why does this method fail for Ton 345? A clue comes from considering the agreement of 
SD SS1228's gravit a tiona l redshi fts as measured from the two different methods. As modeled 
by lGansicke et al\ (l2006l . l2008a] ). the gas disk orbiting SDSS1228 has negligible eccentricity 
while for Ton 345 the disk eccentricity is large (e~0.2— 0.4). In the case of a large disk 
eccentricity it is understandable that measuring the peak midpoint velocity as described 
in Section 13.21 would provide an inaccurate estimate of the systemic velocity. As such, we 
instead take the discrepancy between the gravitational redshifts as computed through the 
two different methods for Ton 345 as additiona l evidence for the large eccentricity of the 



gaseous disk (e~0.2-0.4, iGansicke et a/.l l2008a[ ). It is noted that these high eccentricity 



values in Ton 345's disk persist over long time-spans. Ton 34 5's disk is like l y pumped up to 
(and maintained at) high eccentricity by a surviving planet ( JLubowlllQQll : iKley &: Dirksen 



20061 ). 



Gas Disk Models 



When observed at high spectr al resolution the relatively c l ean-cu t morphologies of the 
gaseous emission lines presented in iGansicke et al\ ( 120061 . 120071 . l2008al ) resolve into complex 
structures. A common feature of the HIRES observed Ca II IRT emission lines for all three 
white dwarfs is a pro file having the typical peak and wing structure (e.g., see Figure 1 in 
Home &: Marsh! Il986l ) on one side of the emission (red side of the double peaked Ca II 
emission features for SDSS1228 and Ton 345 and the blue side for SDSS1043 in Figure [3]) 
while having a peak with a sharp cut-off in emission (i.e., no wing) on the other side (blue 
side of the double peaked Ca II emission features for SDSS1228 and Ton 345 and the red side 
for SDSS1043 in Figure |3]) of the emission. One could imagine that such an emission-line 
morphology would st ill be well reproduced by rnodehng the disks with families of elliptical 
rings as was done by IGansicke et al\ ( l2006l . 120071 l2008al ) . 



Werner et al\ (|2009[ ) present a physical model to explain the gas emission lines observed 



in SDSS1228's spectrum; their model relies on an active disk that is heated by viscous 
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dissipation of disk energy. We use the HIRES spectra to test the predictions of this model. 
In particular, we compare the HIRES spe ctra covering t he Ca II H, K, and IRT tra nsitions 
(Figure E]) to the predictions in Figure 6 of IWerner et al\ (120091 ). IWerner et al\ ( 120091 ) predict 
that the summed Ca II H & K emission complex should have flux that is a factor of ~5 
stronger than the summed Ca II IRT flux. The Ca II H, K, and IRT emission lines d etected 
in the spectra of SDSS1228 can be used as a robust test of the IWerner et al\ (120091 ) model. 
For each emission line the continuum flux is subtracted and the line flux summed. Extraction 
of the Ca II H emission line flux is not attempted as this complex is contaminated by the 
He absorption line. We instead assume that the emission line flux is the same for both of 
the Ca II H & K transitions and multiply the measured Ca II K emission line flux by 2 as a 
proxy for the sum of the emission line flux for both transitions. Examination of SDSS1228's 
emission line flux for the Ca II H, K, and IRT transitions indicates that the summed Ca II 
H & K emission line flux is ~3.5 times weaker than the summed Ca II IRT emission line flux 
(see Table [9]) . Thus, the observational result differs by a factor of ~17 from that expected in 
the IWerner et al\ ( l2009l ) model. It is noted that the total flux in line emission is ~3xl0~^ 
times the bolometric flux emitted by SDSS1228. 

Assuming the line emission is opt ically thick (as is de termined in the modeling pre- 
sented in the supplemental material to iGansicke et a/.ll2006[ ). we can estimate the gas disk 
temperature from our observations by comparing the emission line flux for Ca II transitions 
in different parts of the optical spectrum. For SDSS1228, the only white dwarf with detected 
Ca II H+K emission, we estimate a gas disk temperature of ~5000 K. Gas disk temperature 
upper limits for Ton 345 and SDSS1043 are reported in Table [Hi 



6.1. Z II Region Model 



The IWerner et al\ ( 120091 ) model does not provide a good match to the HIRES data. As 
such, we briefly describe an alternative model for the disk heating mechanism. We envision 
a model where the gas is photoionized by ultraviolet photons from the white dwarf. The 
idea w ould be akin to an H II region, but since this is a gas that is extremely deflcient in H 
or He ( JGansicke et a/.ll2006l 120071 l2008al Jbl. and the results preseneted herein), it would be 
more aptly described as a Z II region. To determine the steady state temperature of a Z II 
region one must consider three processes: 



(i) injection of energy into the gas from photoionized electrons; 

(ii) energy loss from the gas when electrons recombine with ions, a process for Z II regions 
that would be negligible when compared to the last process: 
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(iii) energy loss from the gas due to optically thick line emission. We assume that the gas 
density is sufficiently high (>>10^ cm~'^) that forbidden line cooling is unimportant 
and that gas instead cools through emission lines like the Ca II IRT. 

We expect that the gas disk scale height will be significantly larger than the dust disk 
scale height (see Section r6.1.ip. As a result, collisional energy transfer from gaseous atomic 



species to dust particles (e.g., iGoldreich fc Kwanlll974J ) will be negligible compared to the 



above three energy balance points. Hence, it is expected that the gas and dust temperatures 
are decoupled. 

Here we pro vide details for a highly idealized Z II region model. We consider a flat 



passive disk (see I Jural l2003bl ) whose thickness is a small fraction of the radius (i?*) of the 
illuminating star. The star is assumed to have effective temperature T* and photospheric 
emission well described by a blackbody By(T^). At radial distance D from the star, we 
compute the effective disk-heating flux of stellar photons above energy threshold hvi with 
frequencies between vi and V2i F~{i'i, 1/2) in the limit that D » R^,, as: 

If there is no gas in the system, then all the incident stellar flux heats the dust so that z/i 
= and 1^2 = 00. In a system with both dust and gas where, for simplicity, it is assumed 
that all gas atoms have the same average ionization potential huj, then the grains absorb all 
the photons with u < uj while the gas absorbs all the photons with v > vj. Therefore, the 
flux heating the dust, -F^g^, is computed from Equation (1) for z/i = and U2 = ui. Since 
we consider environments where hvi » kT^,, most stellar photons are absorbed by grains, 
and the expected dust temperature is nearly the same as if no gas was present. For the 
flux heating the gas, F~^g, we use Equation (1) with z/i = z// and z/2 = 00. Using the Wien 
approximation for the Planck curve we find: 

P,as-^{-R] -tH- / (x' -x/x2)e "dx (2) 



3 V D / c 



XI 



where: 



"' = It: <=" 



From expressions (1) and (2), we find: 
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The value of /ii// is determined from the metalhc constituents in a gaseous disk. Hence, 
knowledge of the gas disk composition is imperative in calculating this aver age energy of a 



photo e lectron and t he distribution of its energy into the gaseous material. iGansicke et al. 



(I2OO6I . 120071 . |2008aj) report emission lines from Ca in all three white dwarfs and Fe in 
SDSS1228 and Ton 345, absorption lines of Mg in all three white dwarfs, and Ca and Si 
absorption in Ton 345. Under the assumption that hea vy elements in the photospheres of 



these white dwarfs were accreted from their disks (see e.g.. lZuckerman et a/.ll2007l : lKlein et al. 



2OIOI ). we can claim that each gas disk holds at least the elements Mg, Fe, Si, and Ca. With 
this mix of elements we calculate an average ionizing potential of ~8 eV. Better knowledge 
of the gas disk composition, and ultimately the energetics of the gas disk, is most likely to 
come from deep, high-spectral resolution spectroscopic observations of these white dwarfs in 
the ultraviolet. For the case of a star with effective temperature T^ = 20,000 K and huj = 
8 eV, we compute F-jFg^^ ^ 5%. 

From the measurements reported in Table [9] for SDSS1228, Ton 345, and SDSS1043 
we see that the line fluxes received at Earth summed over the calcium triplet emission are 
roughly 10 x 10~^^, 4 x 10~^^, and 2 x 10~^^ erg cm~^ s~^, respectively. For these same 
three stars, the fluxes from the dust can be very roughly approximated as uF^, evaluated at 
3.6 fim, or 2 X 10"^'^, 2 x 10~^^, and 3 x 10~^^ erg cm~^ s~^. Therefore, for each of the three 
white dwarfs, the ratio of the flux received in the calcium triplet line emission compared to 
the flux received from the dust is ~5%, as predicted. We take this agreement of the data 
with the model as support of our description of the systems as Z II regions. 

The flux re-radiated by the gas in the disk, F^^^ is: 

F+, = Y^n B^iTgg,) Au (5) 

In this expression, the sum is performed over all important cooling lines, here assumed to 
be the three members of the calcium triplet. Although the disk is orbiting rapidly, we 
assume that each portion of the gas is in vertical hydrostatic equilibrium and that the local 
microscopic gas motions are mainly thermal. Therefore: 
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where vl is the line frequency, in{Ca) is the mean atomic weight of calcium, and Pune is a 
coefficient of order unity. We compute the gas temperature by finding the value of Tgas such 
that: 



^gasi-'-gas) — ^ gas- I ' J 

With considerable simplifications, we compute an analytic solution to Equation (6). We 
assign the same frequency to all three members of the calcium triplet and take Pune = 4. 
Combining the previous equations we find: 



^"' V m{Ca) J c3 



For convenience, we define A such that: 

where hu^ = kT^. Combining the above equations, we find: 

TgUD) = tJ^ - ^InA] . (10) 

Gas temperatures are calculated as a function of distance from the star using the analytic 
solution in Equation (10) and are plotted for specific cases in Figure |9l We assume in these 
calculations that m{Ca) = 40 atomic mass units. For D << 20 R^,, the simple model fails 
because there is insufficient cooling. However, for D >> R^ and Teff = 20,000 K, we see 
that the gas temperature ranges from 6000 K to 3000 K, depending upon the exact distance 
from the star. Such gas disk temperatures are roughly in agreement with those estimated 
from observations of the three gas disk-hosting stars. 

6.1.1. Decoupled Gas and Dust Disks 

In our scenario for tidally-disrupted disks orbiting white dwarfs, the gas temperature 
and the dust temperature are distinct because the gas and dust are spatially separated. The 
dust is assumed to be confined to a fiat disk somewhat analogous to Saturn's ring. Depending 
upon the mass of the tidally-disrupted rocky object and the extent of the radial zone where 
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the material is located, the vertical thickness of the dust may be only 1 cm. In contrast, the 
gas is pictured to be in vertical hydrostatic equilibrium and therefore: 

Pgas = P(0)e-^'/^' (11) 

where p(0) is the midplane density. For gas orbiting at radial distance D from the star of 
mass M^,: 

With D = 10^° cm, Tg^s ~ 6,000 K, M, = 0.6 Mq, and adopting /x = 2.4 x 10-^3 g as 
representative of an ionized gas of heavy atoms, we estimate H ^ 3 x 10^ cm. Thus, the 
bulk of the gas is very far from the dust and their temperatures may be different. 



7. Discussion 

From the imaging and spectroscopic data we have estimated inner and outer radii for 
both the gas and dust disks using as few assumptions as possible to connect the dust and 
gas disk results. Now, the two disk parameter sets are combined to understand how the gas 
and dust are related, and what that can tell one about their origin and evolution. 

Dust disk inner radii are close to what is expected to correspond to the sublimation 



tempe rature for silicate dust particles. Spitzer IRS observations presented in IJura et al. 



( 120091 ) suggest white dwarfs with circumstellar disks have orbiting dust composed of olivine. 



The limited wave length coverage o f these spectra preclude more complete compositional 



analyses (but see iReach et a/.ll2009l ). None of the gas disk- hosting white dwarfs discussed 
herein have existing mid-infrared spectroscopic observations, hence there is some uncertainty 
in what dust species orbit these three white dwarf stars and at what temperature these 
materials will sublimate. The dust outer disk radius for SDSS1228 is consistent with the 
expectation for its Roche radius and in light of similar results for several other dust disk- 



hosting white dwarfs (e.g., iFarihi et a/.ll2009l ). we assume the same is true for Ton 345 and 
SDSS1043. 

Gas disk outer radii are spatially coincident (to within measured and modeled uncertain- 
ties) with dust disk outer radii (Figure [8]). The gas disk inner radius for SDSS1228 and Ton 
345 appears to be farther from the host star than the dust disk inner radius. For SDSS1043 
the gas disk is closer to the host star, although the uncertainty range for its gas disk inner 



-18- 



radius is large (Table |8]). Slightly different dislc inclinations (that are still within the range 
of viable model parameters, see Table [7]) could push the dust and gas disk inner radii into 
better agreement. Lower disk inclination values than those quoted in Table [7] would result 
in gas disk inner radii closer to the host white dwarf star and dust disk inner radii farther 
away from the host white dwarf star and vice-versa. 

Similar outer radii of the dusty and gaseo us di s ks for the three w hite dwarfs can be 
understood through the formalism developed in I Jural (120081 ). IJural (120081 ) suggests that pre- 
existing dusty disk material can sputter entering rocky object material from a white dwarf's 
remnant planetary system. The pre-existing material would have been some rocky material 
that was tidally shredded after falling within the white dwarf's tidal radius. New infalling 
material would tidally shred and then sputter out onto the pre-exisiting disk, populating the 
disk with gaseous material. Steady state flows of rocky material could maintain the dusty 
and gaseous disks that are observed around each of the three white dwarfs. 

The above discussion of disk inner radii suggests that the gas and dust disks are likely 
to be entirely spatially coincident from inner to outer disk. However, the best fit model 
parameters (Table [7]) for SDSS1228's and Ton 345's disk systems would have to be stretched 
to the extremes of their viable ranges to be in agreement. Here we discuss another possible 
way that the dust disks around SDSS1228 and Ton 345 could reside farther from the host 
stars and be more in agreement with the gas disk inner radii without changing disk inclination 
angles. Flat dust disk m odels predict the temperatures shown in Table [71 As suggested by 
Brinkworth et al\ (120091 ) . if a warp were present in a dust disk then its grains could be cooler 
and located farther from their host white dwarf star while still reproducing the observed 
excess infrared emissio n. Warps for the dust di sk-hosting white dwarfs GD 362 and GD 56 
have been suggested by I Jura et al\ (l2007bll2009l ). A sligh t warp in the dus t disk of SDSS1228 
or Ton 345 could inflate the inner radius by a factor >1 (jjura et a/.ll2009[ ). Such an inflation 
could bring SDSS1228's and Ton 345's inner dust disk radii into agreement with their inner 
gas disk radii. For the case of SDSS1043, the combined uncertainties for its gas and dust disk 
inner radii are compatible with both terminating at the same inner radius. The asymmetric 
gas emission lines observed (Figure [3]) suggest that the white dwarfs' gas disks are not likely 
to be fiat, lending support to a warped dust disk interpretation. 

Increased disk viscosity at the gas and dust disk inner termination radius can provide 
an explanation for similar gas and dust disk inner radii as well as the lack of emission 
corresponding to material extending inwards to the white dwarf stars. In the region of the 
disk where dust grains are not yet sublimated, the net viscosity of the disk is small and 
disk material does not drift inwards rapidly. The inner disk is heated to the sublimation 
temperature where grains begin to vaporize resulting in inner disk material that is all gaseous 
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and ionized. Higher densities of hot gas leads to enhanced disk viscosity which causes 
material to drift inwards towards the host white dwarf star on short timescales. Rapid inner 
disk depletion will result in small values of inner disk surface density, consistent with the lack 
of any emission having velocities corresponding to radii closer to the white dwarf than the 
observed inner disk radii. Thus, the sublimation of dust grains and resulting enhanced gas 
density can potentially explain coincident gas and dust disk inner radii around these white 
dwarfs. The eventual fate of the inner disk material will be to accrete onto the host white 
dwar f star, polluting its atmosphere and providing a unique insight into its bulk composition 
(e.g., 



Zuckerman et a/.l 120071 : iKlein et a/.ll2010l : iDufour et al.\\2Qm 



8. Conclusion 

We have obtained a suite of ground-based optical and near-infrared measurements of 
three gas-disk hosting white dwarfs. Unambiguous infrared excess emission is identified 
towards Ton 345, confirming that this white dwarf hosts a dusty disk in addition to its 
gaseous disk. Characterization of these three gaseous and dusty disk systems indicates that 
the gas and dust disks are spatially coincident at their outer radii and likely their inner 
radii. Disk parameterization results are consistent with a scenario where both dust and gas 
disks have their origin in the dissolution of rocky objects from the white dwarfs' remnant 
planetary systems. 

Detection of the Ca II H & K and Ca II IRT lines in emission towards SDSS1228 enables 
the development of a new model for the gas disk heating mechanism. This model relies on 
photoionization of metallic atoms in a metal-dominiated region around the hot white dwarfs 
and cooling of the gaseous material through optically thick emission lines. 
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Fig. 1. — Top panels: CFHT WIRCam images of Ton 345 in the J-, H-, and Kg-bands (from 
left to right) with a hnear stretch. The image orientation and scale is the same for all three 
bands. Note the fainter source located ~2.0" to the NE of Ton 345. This faint background 
source does not contribute a significant amount of flux to the CFHT measured magnitudes of 
Ton 345. Bottom panels: Gemini North NIRI images of Ton 345 in the J-, H-, and Kg-bands 
(from left to right) presented on a linear stretch. Images were smoothed over with a 7-pixel 
boxcar to enhance the faint background object contrast. The image orientation and scale is 
the same for all three bands. The background source, presumably a galaxy, is now clearly 
detected in the H- and Kg-bands. 
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Fig. 2. — Gemini North NIRI image of Ton 345 in the L'-band presented on a hnear stretch 
after being smoothed over with a 7-pixel boxcar. The fluctuations in the background around 
Ton 345 are representative of the la noise level. The background galaxy is not detected, it 
is likely too diffuse for the sensitivity threshold of this imaging set. 
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Fig. 3. — Gas emission lines detected in HIRES optical spectra for SDSS1228, Ton 345, and SDSS1043. Top 
Panel: Ca II infrared triplet region. Combined spectra from the 2008 November run are plotted. The lack of 
data around ^8600 A is due to the gap between red orders for the HIRES setup used. Middle Panel: Fe II 
A5169 region. Combined spectra from the 2008 November run are plotted. Emission lines are seen only in 
SDSS1228 and Ton 345. SDSS1228 appears to have three emission peaks where one may be due to emission 
from one of the Mg I b triplet transitions. Bottom Panel: Ca II K region. Plotted spectra are from the 2007 
May run for SDSS1228 and SDSS1043, and from the 2008 February run for Ton 345. Emission lines are seen 
only in SDSS1228 (Ca II H emission for SDSS1228 is detected but not shown). The weak, blueshifted Ca II 
absorption components are likely interstellar in nature based on comp arison to known int erstellar Ca II K 
line strengths and velocities in the direction towards each white dwarf ( Albert et a/.lll993l) . Wavelengths in 
this figure are corrected to the heliocentric reference frame and are presented in air. 
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Fig. 4. — Multiple epochs of gas emission lines detected in the HIRES spectra. Top Panel: Fe II 
A5169 emission lines for SDSS1228. The 2007 May data are from a single exposure. The 2008 Nov 
data are three combined exposures. The S/N of each are roughly equivalent (see Tabled]). The 
relative strength of the three emission peaks seems to vary between the two epochs, although the 
S/N of each epoch prevents any conclusive analysis. Middle Panel: Ca II infrared triplet emission 
lines for Ton 345. All measurable parameters for the two epochs agree to within their respective 
errors. Bottom Panel: Fe II A5169 emission lines for Ton 345. All measurable parameters for the 
two epochs agree to within their respective errors. Wavelengths in this figure are corrected to the 
heliocentric reference frame and are presented in air. 
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Fig. 5. — SED for the DBZ white dwarf Ton 345. Overplotted is a model for the white 
dwarf photospheric emission (a blackbody with Tefr~18,600 K represented by the dotted 
line) and an orbiting fiat, passive, opaque dust disk (represented by the dashed line). Broad- 
band fluxes can be found in Table [1] and the model parameters can be found in Table [3 
Horizontal bars correspond to filter bandpasses. Optical fluxes were obtained from SDSS 
while ultraviolet fluxes were taken from the GALEX AIS and MIS catalogs. Kg and L' are 
in excess of what one would expect from the photosphere of the white dwarf alone. 
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Fig. 6. — SED for the DAZ white dwarf SDSS1228. Overplotted are models for an at- 
mosphere of a Teff=22,020 K, log(5f)=8.24, DA white dwarf (dotted line; D. Koester, 2009 
private communication) with an orbiting flat, passive, opaque dust disk (dashed hne). Input 
infrared fluxes are taken from iBrinkworth et al\ (120091 ). Three different disk models are dis- 
played to illustrate a range of viable model parameters. The black model curve (the middle 
of the three) is the best fit and the model whose parameters we adopt for SDSS1228's dusty 
disk (see Table [7j). Optical fluxes were obtained from SDSS while ultraviolet fluxes were 
taken from the GALEX AIS catalog. One possible explanation for the r elatively high flux 
at 16 /im is that the source possesses a strong silicate emission feature ( jBrinkworth et al. 
20091 ). thus it is not included in the continuum fit. 
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Fig. 7. — SED for the DAZ white dwarf SDSS 1043. Overplotted is a model atmosphere for 
a Teff=18,330 K, log(5')=8.09, DA white dwarf (dotted hne; D. Koester, 2009 private com- 
munication) with an orbitin g fiat, passive, opaqu e dust disk (dashed line). The model UV 
fluxes are reddened using the lCardelli et all ( 1l989l ) extinction curve assuming E(B— V)=0.03. 
Fluxes can be found in Table |2] and the model parameters are presented in Table [71 Data 
points longward of 2 fim are IRAC measurements reported in C. Brinkworth et al. (2010, 
in preparation). Optical fluxes were obtained from SDSS while ultraviolet fluxes were ob- 
tained from the GALEX AIS catalog. When plotted with the Spitzer IRAC data, the Kg 
measurement suggests an excess beginning at ~2 /im. 
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Fig. 8. — Sketch illustrating the results presented in Table [HI The dimensions of each white 
dwarf's dusty and gaseous disks are plotted in units of each white dwarf's radius. True 
blue colored regions correspond to gas disks while cardinal red colored regions correspond to 
dust disks. The vertical scale heights of the gas disks are expected to be larger than those 
of the dust disks (Section 16.1. ip . Gold, horizontally hatched regions correspond to the la 
uncertainty for the dust disk inner and outer radii as reported in Table [HI Gold, vertically 
hatched regions correspond to the la uncertainty for the gas disk inner radii as reported in 
Table E 
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Fig. 9. — Disk temperature as a function of radial separation from its host white dwarf 
star in units of stellar radii. The solid lines are for the gas disk temperature as calculated 
using the analytical Z II region model outlined in Section 16.11 Dashed lines are dust disk 
temperatures as described in Section HI The red, green, and blue curves are for stellar 
effective temperatures of 20,000 K, 15,000 K, and 10,000 K respectively. We note that the 
gas disk model predictions fail for D/R^, < 20 (see Section [6?T|) . 
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Table 1. Broad-band Fluxes for Ton 345 



Band 


A 
nm 




mag 






?lux Density 
(mJy) 










NIK 










WIRCam^ 


NIRji^ 


Average 


Fobs 


F ^ F 

-L * ^ excess 


L' 


3761 


- 


15.2±0.3 


- 


0.21±28% 


0.075 0.135±0.059 


Ks 


2157 


15.95±0.03 


16.04±0.03 


15.99±0.02 


0.27±2% 


0.19 0.076±0.005 


H 


1654 


16.18±0.05 


16.33±0.02 


16.26±0.07 


0.33±7% 


0.29 0.036±0.023 


J 


1215 


16.36±0.03 


16.41±0.03 


16.39±0.02 


0.45±2% 


0.45 0.00 


SDSS DR7- Optical 


z 


893.1 




16.615±0.018 






0.82±1.8% 


i 


748.1 




16.348±0.015 






1.05±1.5% 


r 


616.5 




16.077±0.014 






1.35±1.4% 


9 


468.6 




15.728±0.035 






1.86±3.5% 


u 


355.1 




15.559±0.015 






2.25±1.5% 








GALEX^ 


— Ultraviolet 






NUV 


227.1 




15.50±0.03 






2.28±3% 


FUV 


152.8 




15.69±0.05 






1.92±5% 



'^F^, is the predicted photospheric flux density in the given bandpass assuming Ton 345 
is a blackbody with Teg of 18,600 K and that the J-band monochromatic flux is entirely 
photospheric in nature. 

WIRCam and NIRI magnitudes are in the Johnson syster n. GALEX measureme nts are 
in AB magnitudes. GALEX uncertainties are as suggested in lMorrissev et all (|2007l ). 
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Table 2. Broad-band Fluxes for SDSS1043 



Band 


A 


mag 


-t* obs 




nm 




(mJy) 






WIRCam/Gemini^'''- 


NIR 


Ks 


2157 


17.69±0.05 


0.056±5% 


J 


1210 


17.72±0.14 


0.133±13% 


SDSS DR7 - Optical 


z 


893.1 


17.942±0.036 


0.24±3.6% 


i 


748.1 


17.701±0.017 


0.30±1.7% 


r 


616.5 


17.433±0.014 


0.39±1.4% 


9 


468.6 


17.123±0.022 


0.51±2.2% 


u 


355.1 


17.358±0.022 


0.43±2.2% 


GALEX^ - Ultraviolet 


NUV 


227.1 


17.07±0.03 


0.54±3% 


FUV 


152.8 


16.86±0.05 


0.66±5% 



'^Kg-band data are from CFHT WIRCam 
observations while J-band data are from Lick 
Gemini observations (see Section |2]). 

^WIRCam and Gemini magnitudes are in 
the Johnson system. GALEX measurements 
are in AB magnitud es. GALEX unc e rtaint ies 



are as suggested in iMorrissey et al\ (120071 ) 
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Table 3. NIRI Fluxes For Background Galaxy Companion to Ton 345 



Band mag Flux Density 



J 


>21.0 


<6.5 


H 


19.15±0.24 


22.9±22% 


Ks 


18.47±0.15 


27.3±14% 



Note. — The background galaxy 
at 1.8" separation from Ton 345 ap- 
pears to be a J-band drop-out galaxy. 
Comparing to galaxy models (e.g 
see Figure 11 of 



Sanders et al. 2007 



and references therein) suggests that 
this is a Balmer-break Spiral or Cool 
ULIRG galaxy at a redshift 2~2. 



Table 4. HIRES Observations Summary 



UT Date 


Setup 


Coverage 


Integration Time (sec) 


S/N^ 


SDSS1228 










05 May 2007 


UV Collimator 


3120-5950 A 


3000 


22 


14 Nov 2008 


Red Collimator 


4500-9000 A 


1800 


12 


15 Nov 2008 


Red Collimator 


4500-9000 A 


1600 


13 


16 Nov 2008 


Red Collimator 


4500-9000 A 


1900 


12 


Ton 345 










13 Feb 2008 


UV Collimator 


3130-5960 A 


2x1800 


38^ 


14 Feb 2008 


UV Collimator 


3130-5960 A 


2x1500 


38^ 


26 Feb 2008 


Red Collimator 


4600-9150 A 


3600 


42 


14 Nov 2008 


Red Collimator 


4500-9000 A 


2x1800 


45^ 


16 Nov 2008 


Red Collimator 


4500-9000 A 


1800 


32 


SDSS1043 










05 May 2007 


UV Collimator 


3120-5950 A 


3000 


13 


15 Nov 2008 


Red Collimator 


4500-9000 A 


2x2000 


18^ 


16 Nov 2008 


Red Collimator 


4500-9000 A 


1500 & 1700 


14b 



00 



Note. — A resolving power of ~40,000 was measured for all epochs from the 
FWHM of single arclines in comparison spectra. 

'^S/N measured at 5750 A. 



'S/N for combined exposures. 



00 



Table 5. Gas-disk Emission Line Measurements 



Transition 


UT Date'' 


Equivalent Width 


Peak Separation'' 


Full Width<^ 


Peak Midpoint Velocity'' 


VmaxSini'' 






(A) 


(km s-i) 


(kms~i) 


(kms-i) 


(km s-i) 


SDSS 1228 














CallK 


05 May 2007 


0.6±0.1 


707±20 


1014±44 


-35±10 


455±23/559±38 


CallH 


05 May 2007 


_f 


731±18 


1020±53 


-22±9 


489±38/531±38 


Fe II A5169'^ 


05 May 2007 
14,15,16 Nov 2008 


1.3±0.1 
1.3±0.1 


_ 


_ 


_ 


_ 


Ca II A8498 


14,15,16 Nov 2008 


13.8±1.0 


656±9 


938±16 


-22±4 


379±7/560±14 


Ca II A8542 


14,15,16 Nov 2008 


19.3±1.4 


655±9 


944±9 


-17±4 


393±7/551±14 


Ca II A8662 


14,15,16 Nov 2008 


19.8±1.7<= 


659±5 


933±35 


-18±3 


382±7/548±35 


Ton 345 














Fe II A5169 


13,14.26 Feb 2008=* 


1.1±0.2 


838±14 


1349±87 


-79±14 


- 


14,16 Nov 2008 


1.1±0.2 


869±18 


1305±96 


-69±18 


- 


Ca II A8498 


26 Feb 2008 
14,16 Nov 2008 


5.2±0.4 
4.9±0.4 


842±20 
827±18 


1305±49 
1252±73 


-90±11 
-104±8 


671±35/634±35 
619±18/634±71 


Ca II A8542 


26 Feb 2008 
14,16 Nov 2008 


4.1±0.3 
3.5±0.3 


857±16 
891±13 


1281±63 
1211±56 


-98±8 
-89±6 


706±35/575±53 
618±18/593±53 


Ca II A8662 


26 Feb 2008 
14,16 Nov 2008 


3.4±0.2 
3.4±0.5'= 


876±17 
882±25 


1298±75 
1263±71 


-99±9 

-89±12 


663±28/635±69 
611±17/652±69 


SDSS1043 














Ca II A8498 


15,16 Nov 2008 


7.0±0.5 


757±38 


1482±50 


39±19 


786±35/686±35 


Ca II A8542 


15,16 Nov 2008 


7.0±0.6 


778±49 


1439±50 


52±25 


811±35/628±35 


Ca II A8662 


15,16 Nov 2008 


8.1±1.2'= 


806±48 


1402±109 


42±24 


802±104/600±35 



00 
00 



''Spectra have been combined when multiple days are listed. 

''See Sections 13.21 and [5l for description of these quantities. The two different values reported for 
velocity gas seen in the blue and red wings of the double-peaked emission features, respectively. 



xsim correspond to the maximum 



'^FuU velocity width of emission feature, from continuum blueward of the blue emission peak to continuum redward of the red emission 
peak. 

''To within the noise of each individual exposure the 13,14 Feb 2008 and 26 Feb 2008 Fe II A5169 emission features for Ton 345 are 
identical; to improve the S/N we combine all three spectra into a single 2008 Feb epoch. 

'^Measurement may be contaminated by second order blue light leak, see Section 12.41 

Photospheric Balmer He absorption in SDSS1228 prevents an accurate continuum estimation for the Ca II H emission complex. 

SThe Fe II A5169 emission in SDSS1228 shows peculiar features (see Figure |4]|, the origin of which are uncertain. As such we do not 
include these measurements in this table. 



00 

CO 



Table 6. Kinematics of Gas-disk White Dwarfs 



WD Name RVobs'' Mwd Rwd Vgrav*" 
(kms-i) (Mo) (Rq) (kms-i) 



Distance pmRA pmDE U V W 

(pc) (mas yr~^) (mas yr""'^) (km s~^) (km s~^) (km s~^) 



SDSS1228 


+37±4 


0.77 


0.011 


+44 


125+13 


-44.0+2.0 


-18.0+1.0 


-18 


-20 


Ton 345 


+36±4 


0.70 


0.010 


+44 


85+10 


-10.0+1.0 


-62.0+2.0 


+9 


-21 


SDSS1043 


+39±4 


0.67 


0.012 


+35 


185+19 


+38.0+6.0 


-46.0+4.0 


+45 


-27 



-12 
-14 
f6 



''RVobs is the observed radial velocity as measured from the absorption lines of each white dwarf. It has not been corrected 
for each white dwarf's gravitational redshift. 

'^Gravitational redshift as computed from the model-determined mass and radius reported for each white dwarf in columns 3 
and 4. 



o 



Note. — See Section 13.21 for discussion of input and output parameters. For the UVW space motions positive U is towards 
the Galactic center, positive V is in the direction of Galactic rotation, and positive W is toward the north Galactic pole. 
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Table 7. Dust Disk Models 





SDSS1228 Ton 345 


SDSS1043 




Stellar Parameters 


Teff (K) 

Rwo/Dist 


22,020 18,600 
2.0 X 10-12 2.6 X 10-12 


18,330 

1.5 X 10-12 


Disk Parameters'^ 



Tinner (K) 1100 1300±50 1350 1500 1200 

Touter (K) 700 500±70 350 1000 1100 

i (°) 60 73±3 76 66 60 

^The three different parameter sets for SDSS1228 correspond 
to the three curves plotted in Figure O The parameters are for 
the red, black, and blue curves respectively. We adopt the middle 
parameter set (having Tinner of 1300 K) for SDSS1228's dust disk 
properties (corresponding to the black curve in Figure [6]). From 
the three curves we estimate the la uncertainties on the model 
parameters (see Section H]) which are reported with SDSS1228's 
best fit parameters. 

Note. — The effective temperatures for the three white dwarfs 



come from iGansicke et all (120061 120071 l2008al ). The ratio of the 



radius of the white dwarf to the distance of the white dwarf from 
Earth, Rwo/Dist; temperature at the inner edge of the disk. 
Tinner! cLud inclination angle, i, come from the model fit to the 
U V, optical, and infrared photometry. See Section H] for a discus- 
sion of the temperature at the outer edge of the disk (Toutor)- 
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Table 8. Disk Dimensions 







Gas Disks 






Dust Disks 




Name 


^max 


dinner, gas 


^outer,gas 


^inner,dust ^outer,dust 


J- outer 




(km s-i) 


(Rwd) 


(Rwd) 


(Rwd) 


(Rwd) 


(K) 


SDSS1228 


575±17 


40±3 


108 


26±2 


93±17 


500 


Ton 345 


709±20 


27±2 


~100 


17±2 


100±18 


400 


SDSS1043 


923±52 


12±2 


81 


23±2 


80±15 


470 



Note. — All quoted uncertainties are la (see below). Gas disks: Vmax comes 
from the average of the blue or red wing values (whichever is higher) listed in 
the "vinaxsini" column of Table [5l This average value is then divided by the sine 
of the inclination angle reported in Table El l^inner,gas is derived directly from 
the value reported in the "vmax" column. The quoted uncertainties for Rinner,gas 
also include the statistical uncertainties for the white dwarf masses and radii as 



quoted in iGansicke et al\ (J2006l . 120071 . l2008al ) . There may be additional syst e matic 
unce r taintie s that are not included. Router,gas is taken from IGansicke et al\ (J2006l . 
20071 . l2008aH bl) . Dust disks: See Section H] for a description of how Rjnner,d«st and 
Router,dust are derived. Touter corresponds to the outer dust radius shown. The 
uncertainties for Rinner,dust and Router,dust are estimated from the range of possible 
disk parameters modeled for the SDSS1228 dat a set and the suggested uncertaint y 
for each star's effective temperature as quoted in lGansicke et al\ (J2006l . 120071 . l2008al ) . 



Table 9. Ca II Emission Line Flux 



Name Ca II K Ca II A8498 Ca II A8542 Ca II A8662 Ratio Gas Disk 

Firt/Fhk Temperature (K) 



Total Line Flux*^ (10 ^^ ergs cm ^ s ^) 



SDSS1228 


1.4 


2.6 


3.6 


3.7 


-3.5 


-5000 


Ton 345 


<2.b'° 


1.8 


1.3 


1.2 


>0.9'^ 


<6500 


SDSS1043 


<0.6'^ 


0.7 


0.7 


0.8 


>2.0'' 


<6000 



^These values are computed by multiplying emission line EW measurements reported in Table [5] 
by the stellar continuum flux at the emission line location. 

''These upper limits are calculated assuming we could not detect Ca II H & K emission features 
weaker than those seen in SDSS1228 (see Figure[3]). 



00 



Note. — See Section [6] for discussion. 



